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Structure function evolution at next-to-leading order and beyond 

Andreas Vogt''* 

'^Instituut-Lorcntz, University of Leiden, P.O. Box 9506, 2300 RA Leiden, The Netherlands 

Results are presented of two studies addressing the scaling violations of deep-inelastic structure functions. 
Factorization-scheme independent fits to all ep and fip data on F2 are performed at next-to- leading order (NLO), 
yielding as{Mz) = 0.114 ± 0.002cxp (-1-0.006 — 0.004)th . In order to reduce the theoretical error dominated by the 
renormalization-scale dependence, the next-higher order (NNLO) needs to be included. For the fiavour non-singlet 
sector, it is shown that available calculations provide sufficient information for this purpose at x > 10~^. 



1. Introduction 

One of the important objectives of studying 
structure functions in deep-inelastic scattering 
(DIS) is a precise determination of the QCD scale 
parameter A (i.e., the strong coupling as) from 
their scaling violations. In this talk we briefly 
present results of two studies [ |l|, g| aiming at 
an improved control and a reduction of the cor- 
responding theoretical uncertainties. 

2. Flavour-singlet evolution in NLO [ |l[| 

The evolution of structure functions is usually 
studied in terms of scale-dependent parton den- 
sities and coefficient functions. In this case the 
predictions of perturbative QCD are affected by 
two unphysical scales: the renormalization scale 
fXr and the mass- factorization scale fif. While 
the former is unavoidable, the latter can be elim- 
inated by recasting the evolution equations in 
terms of observables [ y]. In the flavour-singlet 
sector, this procedure results in 
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with Fb ^ dF2/d\nQ'^ or Fb ^ Fl- The ker- 
nels V are combinations of splitting functions and 
coefficient functions which become prohibitively 
complicated in Bjorken-x space at NLO. Thus 
Eqs. (yj) are most conveniently treated using mod- 
ern complex Mellin-moment techniques [ m . 



We have performed leading-twist NLO fits to 
the FP data of SLAG, BCDMS, NMC, HI, and 
ZEUS. Statistical and systematic errors have been 
added quadratically, the normalization uncertain- 
ties have been taken into account separately. The 
singlet/non-singlet decomposition has been con- 
strained by the F^/F^ data of NMC. The ini- 
tial shapes F2,b{x, Qq) are expressed via standard 
parametrizations for parton densities at ^f = Qo- 
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Figure 1. The dependence of the fit results for the 
energy- momentum sum and for as{Mz) on the 
Q^-cut imposed in addition to W'^ > 10 GeV^. 



In order to establish the kinematic region which 
can be safely used for fits of a^ in the leading- 
twist NLO framework, the lower Q^-cut applied 
to the data has been varied between 3 and 30 
GeV^. When the normalized momentum sum of 
the partons defining the i^2,B initial distributions 
is left free, the fits with Ql^^ < 10 GeV^ prefer val- 
ues significantly different from unity, see Fig. 1. 
Also shown in this figure is the Q^^^-dependence 
of the fitted values for as{Mz), now imposing 
the momentum sum rule. The results for Q^ut — 
7 GeV^ tend to lie above the Q'^^^ > 10 GeV^ 
average of as{Mz) = 0.114 (dashed line). 

In Fig. 2 we display the renormalization scale 
dependence of the as{Mz) central values for the 
safe choice Ql^t = 10 GeV^. The conventional, 
but somewhat ad hoc, prescription of estimat- 
ing the theoretical error by the variation due to 
0.25 < Mr / Q^ < 4 results in 
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Other theoretical uncertainties are considerable 
smaller and can be neglected at this point. The 
uncertainty due to possible higher-twist contribu- 
tions, for instance, can be estimated at about 1% 
via the target-mass effects included in the fits. 
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Figure 2. The dependence of the optimal values 
for as{Mz) on the renormalization scale /i^. 



3. Non-singlet evolution in NNLO [ §] 

The theoretical error in Eq. (0) clearly calls for 
NNLO analyses. The necessary contributions to 
the /3-function [ and the coefficient functions 
[ o are known. However, only partial results 
are available for the 3-loop terms p(^^(x) in the 
splitting- function expansion (a^ = as/47r) 



P = a,pW -I- a^pd) + a3p(2) 



(3) 



For the non-singlet part of F2 considered here 
(NS+), present information comprises the lowest 
five even-integer moments [ M , the full A^? piece 
[ ||] , and the most singular small- a; term [ ^ . 

We have performed a systematic study of the 
constraints imposed on P^^ (x) by these results. 
Four approximations spanning the current uncer- 
tainty range are shown in Fig. 3, together with 
their convolutions with a typical input shape. 
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Figure 3. Representative approximate results 
for the fiavour-number independent part of the 
3-loop non-singlet+ MS splitting function. 
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Figure 4. The first three steps in the expansion 
of the scaling violations of the non-singlet com- 
ponent of F2 for typical input parameters. 

P^g {x) is well determined at a; > 0.15, with 
a total spread of about 15% at a; ~ 0.3. At (non- 
asymptotically) small a: its behaviour is rather 
unconstrained despite the known leading a; — > 
contribution. As the splitting functions enter 
scaling violations always via convolutions 



(P®/)(x)= / dy/yP{x/y)f{y) 
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with smooth initial distributions /(x), the resid- 
ual uncertainties are much reduced for observ- 
ables over the full x-range. In the present case 
they prove to be negligible at a; > 0.02. 

The net effect of the NNLO correction is finally 
illustrated in Fig. 4, where the scale-derivative of 
F|^^ is shown for fir = Q and Nf = 4, using an 
a^-value typical for the fixed-target region. The 
inclusion of this correction into fits is expected to 
lead to a slightly lower central value for as and a 
considerably reduced theoretical uncertainty. 

4. Summary and outlook 

We have analyzed present ep/ fip _F2-data in a 
factorization-scheme independent framework [ m . 
We find that Q^, W^ > 10 GeV^ is a safe region for 



leading- twist NLO fits of ag. Our central value is 
close to that of the standard pre- HERA analysis 
in [ |l^ , but lower than the recent result of [ 11 
using a lower Q^-cut of 2 GeV^. The irreducible 
renormalization-scale uncertainty turns out to be 
larger than expected from [ |T^ . 

We have derived approximate a;-space expres- 
sions for the 3-loop non-singlet splitting functions 
P]^g, including error estimates [ ||]. This ap- 
proach is complementary to, but more flexible 
than, the integer-moment procedures pursued in 



[ 12, 13]. The remaining uncertainties of Pj^g are 



(2) 



small for the evolution at a; > 10~ , thus allow- 
ing for detailed NNLO analyses in this region. An 
extension to the singlet case is in preparation. 
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